Low-frequency spectroscopy of superconducting photonic crystals 
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Transmission, reflection and absorption of electromagnetic radiation and photon dispersion law for 
2D photonic crystals with superconducting elements are studied. The calculation of optical proper- 
ties of photonic crystals is studied by layer-by-layer Korringa-Kohn-Rostoker techniques. The results 
of numerical calculations performed for an array of superconducting cylinders have been understood 
in terms of a simple analytical model. The controlling of optical properties of superconducting 
photonic crystal by temperature and magnetic field is discussed. The variation of superconducting 
component density with temperature leads, particularly, to significant reduction of the transmission 
peaks and even to a nonmonotonous behavior of some absorbance peaks. 
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I. INTRODUCTION 



II. PROBLEM DEFINITION AND METHOD OF 
CALCULATION 



Photonic crystals (PCs) consisting of metallic wires 
are a subject of growing interest during last years and 
are widely studied both theoretically [H@13|3@l and ex- 
perimentally [10] • At low- frequencies, e.g., in the giga- 
and terahertz spectral regions, the real part of the di- 
electric function of metallic component is large and neg- 
ative which allows the formation of the complete two- 
dimensional photonic band gap in these spectral regions. 
Moreover, the electromagnetic waves cannot propagate 
in such a system at frequencies lower than the so-called 
cutoff frequency w c ■ ^ lattice constant of these struc- 
tures lies in the millimeter and submillimcter ranges their 
fabrication is relatively simple. However, the damping of 
electromagnetic waves in metals can suppress many po- 
tentially useful properties of metallic PCs. The possible 
solution of this problem is the replacement of the PC 
metallic component by a superconducting (SC) compo- 
nent. The dielectric function depends on superconduct- 
ing gap A (see, e.g., Rcfs. 8,9] and references therein) 
and SC state can be varied by alteration of external 
parameters, such as temperature and external magnetic 
field which provide a method for a control of the optical 
properties of metallic PCs. The dielectric losses are re- 
duced essentially below superconducting transition in SC 
state, and recent experiments have shown that photonic 
band edges become sharper in SC metals 0. 



The band structure of a superconducting PC was an- 
alyzed in Ref. Q. Here we calculate the transmission, 
reflection, and absorption spectra and analyze control- 
ling of these spectra by temperature. 



We study the photonic crystal formed by a periodic 
array of parallel superconducting cylinders with radius 
R arranged in a square lattice with the lattice constant 
a, see Fig. 1. Electromagnetic waves propagating in this 
structure are described by Maxwell equation 



VxVxE(r)=(-] e(r,u)E(r) 



(1) 



where e(r, u>) is the dielectric function equal to that of the 
superconducting metal inside the cylinders and to unity 
outside them (in vacuum). For simplicity, we consider 
the electro-magnetic response of superconductor cylin- 
ders by using the two-liquid Gorter-Kazimir model and 
describe the electronic system as an admixture of two in- 
dependent carrier liquids, superconducting and normal. 
With the increasing of temperature T or magnetic field 
B, the density of the normal component, n n , increases 
at the expense of the superfluid component, n s , because 
the total density of electrons n tot = n„ + n s is conserved. 
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FIG. f : Schematic representation of a two-dimensional pho- 
tonic crystal with cylinders of the radius R arranged in a 
square lattice with the lattice constant a. 
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At the critical temperature T = T c only normal phase 
survives whereas, at zero temperature, only superfluid 
phase does exist. 

In the Gorter-Kazimir model the dielectric function of 
a superconductor can be presented in the following simple 
form 

e GK (u) = 1 - u\ — + . - . . (2) 

y \U) Z lo(uj + Z7j J 

where ui p is the plasma frequency (4-7m tot e 2 /to) 1 / 2 cor- 
responding to the total electron density n to t, m is the 
electron effective mass, 7 is the phenomenological damp- 
ing coefficient describing the relaxation of normal elec- 
tron component, and a is the fraction of superconducting 
component n s /n to t which varies continuously with tem- 
perature (at magnetic field H =0) from a = 1 at t = to 
a = at T = T c . It can be controlled also by magnetic 
field. It should be noted, that this equation is valid qual- 
itatively for the frequencies us below the superconducting 
gap 2A. 

At the critical temperature T c or critical magnetic field 
H cr the superconducting phase vanishes and the dielec- 
tric function can be described by the conventional Drude 
equation 

egkM = 1 — rrr^i 

uj{u> + 27) 

In the opposite limit, T — > 0, the contribution of the nor- 
mal component is negligible and the dielectric function 
reduces to 

£gk(w) = 1 1 



In numerical calculations we consider high-temperature 
superconductor YBaCuO; its parameters in the 
normal state are as follows: o; p =1.67x 10 15 rad/s, 
7=1. 34x 10 13 rad/s. For YBaCuO the superconducting 
gap can be estimated as 2A ss QksT c /h (ks is the 
Boltzmann constant). For T c =91 K one has A=28 meV 
which is equivalent to 7 THz. A value of 150 /im is 
chosen for the lattice constant a so that frequencies of 
the first few allowed photonic bands lie below A. The 
filling factor f=n(R/a) in our calculations is equal to 
5%. 

The optical spectra are calculated for finite-size two- 
dimensional PC slab with the normal parallel to the sym- 
metry direction (01). In this case the structure can be 
considered as N identical layers, each containing one- 
dimensional chain of cylinders. In what follows we re- 
strict ourselves to the case of normal incidence of TE- 
polarizcd light, which means that the electric field is 
parallel to the cylinder axis z, E = (0,0, E z ) In this 
particular geometry Eq. ([T]) can be rewritten as 

d 2 d 2 /lu\ 2 



E z (x,y) = (3) 




dx 2 dy 

The numerical calculation is performed by a two- 
dimensional bulk and layer-by-layer Korringa-Kohn- 
Rostoker techniques fTol ] . 
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FIG. 2: Transmission (a), reflection (b), and absorption (c) 
spectra of a superconducting photonic crystal calculated for 
three different values of a and the set of other parameters 
given in the text. At a = 1 the absorption is zero and the 
corresponding line A(u>) = coincides with the abscissa. 



III. RESULTS AND DISCUSSION 

Figure 2 presents the transmission, reflection and ab- 
sorption spectra, respectively, T(ui), R(u) and A(u>) = 
1 — R(ui) — T(co), for a 8-layer-thick superconducting PC 
calculated for different values of parameter a. The chosen 
spectral region covers the first allowed band, diffraction 
channels at these frequencies are closed. For the pure SC 
phase, a = 1, the absorption vanishes, R(u>) = 1 — T(w), 
and the transmission peak values reach unity. With de- 
creasing a the transmission and reflection at any fixed uj 
decrease monotonously, with the transmission being more 
sensitive to the variation of a. At the same time the ab- 
sorption maximum can exhibit a non-monotonous behav- 
ior, at least it is the case for the first absorption peak. 
For the chosen lattice constant the parameter 2-7rc/a is 
comparable with the attenuation 7, the role of the elec- 
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tronic relaxation is significant except for very low fraction 
of normal electrons, i.e. except for the region 1 — a <C 1. 
The oscillations in optical spectra can be interpreted in 
terms of the interference effects taking into account the 
modification of the photonic dispersion in the PCs. Such 
spectral behavior is typical for metallic PCs, see, e.g., [If. 
However, for the THz frequency region considered here, 
the values of |£gk(w)| ar< 3 of the order of 10 4 which is 
by almost three orders of magnitude larger than for the 
optical region studied in 0]. 

Below we propose a simple analytical description of 
photonic modes near the edge of the allowed photonic 
band where the photon frequency and two-dimensional 
Bloch wave vector k = (k x , k y ) are related by 



uj = luq + Ck 



(4) 



Here the zero-fc frequency ojq and the coefficient Q have 
imaginary contributions if Im{eQK(w)} ^ 0. Taking into 
account that at low frequencies real part of £gk(<^) is 
negative and its modulus is very large, we can expand u>o 
and C in powers of 1/ y/— £gk- In the first-order approx- 
imation we obtain instead of Eq. (01 



(5) 

where c is the light velocity in vacuum (introduced here 
for convenience), lu^° and c 2 rj are the values of imq and 
C in the limit £gk — » oo. It should be stressed that the 
parameters fl, r\ and C, depend only on the geom- 
etry of the system and arc independent of the cylinder 
dielectric function egc- More precisely, one can use the 
representation 



2?rc 



2nc 



a a 

= 3i , it = s 2 , 7? = s 3 - — , C = s 4 ^ — > 

a a Zirc Zirc 

where the coefficients si to S4 depend only on the ra- 
tio R/a. The fitting of the first-allowed-band dispersion 
curves obtained for different values of £sc give the fol- 
lowing values 



si 



0.415 , s 2 = 0.63 , s 3 = 0.98 , s 4 = 3.45 



Equation (O can be compared to the corresponding ex- 
pansion for a planar waveguide near the cutoff frequency 
see [Tl[. For an ideal planar waveguide of the thickness 
d, one has 



ttc d . , , 2 
~~d + 2^~c (cfc) 



The cutoff frequency and coefficient r\ are equal to nc/d 
and d/2irc, respectively. On the other hand, for the PC 
under study, the frequency is close to ttc/ a (or, equiv- 
alently, si f=a 0.5) and S3 ks 1. This helps the interpreta- 
tion of Eq. ([5]); namely, in the PC the effective waveguide 
"boundaries" are the planes containing cylinders. When 
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FIG. 3: (a) Transmission spectrum T(u>) and (b) photon dis- 
persion k' (lo) near the first transmission peak. Solid, dashed 
and dotted curves are calculated for a = 1, a = 0.5, and 
I esc I — » co- Other parameters including N — 8 are the same 
as in Fig. 2. The vertical lines indicate the frequencies at 
which k'(uj) = ir/(Na), this value of k' is indicated by a hor- 
izontal line. 



the waveguide is not ideal, i.e., the dielectric constant 
is negative, large but finite, the radiation slightly pene- 
trates in its walls which leads to an increase in the effec- 
tive thickness d and a decrease in the cutoff frequency, in 
agreement with the positive sign of f2. 

The relation between the transmission peak position 
and the photon dispersion is analyzed in Fig. 3. Figure 3a 
presents the first transmission peak for the PC structures 
with a = 0.5, a = 1 and esc ~^ o° ■ The first two lines 
are fragments of those shown in Fig. 2a. The third line 
corresponding to esc —* °o (equivalent to the boundary 
condition E z (p) = at the cylinder surfaces) looks like 
a slightly blue-shifted spectrum of the purely SC PC. 
A value of the shift is in agreement with the estima- 
tion Q [— £gk(w§°;7 = O)] -1 / 2 that follows from Eq. ([5]). 
This equation also allows one to understand the evolu- 
tion of the spectrum with decreasing a: for a different 
from unity the photon wave vector acquires an imaginary 
part and the transmission decreases. The similar analysis 
clearly shows that the spectral position of the transmis- 
sion maxima in Fig. 2 are determined by the Fabry-Perot 
interference condition 



Nk'{uj m )a = mir, m = 1, 2 . . . N - 1 , 



(0) 



where k'(u> m ) is the real part of the photon wave vector 
at the frequency ui m . Thus, the peaks in the transmission 
spectrum can be related to the Fabry-Perot interference. 
Extending this one-dimensional interpretation we can ap- 
ply the approximate estimation for the transmission and 
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reflection coefficients 



T n (uj) 



e -1k"(u)aNU _ r 2|2 



Rn(co) = 



M _ r 2 e 2ifc(w)aJV|2 ' 
| r n — e 2ifc(w)a7V-j|2 



(7) 



h _ r 2 e 2ife(w)aAf |2 



derived for a slab of the thickness Na. Here fc" is the 
imaginary part of the wave vector k{u>) calculated for 
the photonic crystal and r is the internal reflection co- 
efficient considered as a fitting parameter while making 
a comparison with the exact calculation. Equation (J6j) 
implies the condition r" = and, for lowest peak with 
m = 1, the fit gives r r* 0.84. 

The two-dimensional nature of the PC makes it diffi- 
cult to derive a stricter analytical equations for the spec- 
tra, but the key features of their dependence on a, in par- 
ticular, the red shift and decrease of the absolute value 
of the transmission are quite well reproduced by Eqs. ([7]) 
with fc(cj) determined from Eq. ([5|). Fig. 4 presents the 
same curves as Fig. 2 but in large scale and only in the 
vicinity of the peak m = 1. One can see that dashed 
curves, calculated by the approximate equations (J7J) lie 
sufficiently close to the exact ones. Expressions (|7|) can be 
simplified in the spectral region near each frequency u) m 
satisfying Eq. ©. Assuming 1 - r « 1 and k"aN -C 1, 
we can reduce Eqs. (|7|) 



,(0)2 



A n (oj) 



or r (0) 

(w - ^,„) 2 + (r^ + r m ) 2 



(8) 



with 



T m = k"{w m )aNA, = (1 - r)A, 



A = 



dk'(u> r , 



(9) 



It follows then that narrow peaks in the transmission and 
absorption as well as narrow dips in the reflection can be 
described in terms of resonant tunneling in the region of 
a "confined" -photon eigenfrequency uj m . Consequently, 
Tm and T m mean the radiative and non-radiative de- 
cay rates of the corresponding "confined" state m. The 
non-radiative decay vanishes for a — 1 and increases 
with decreasing a which leads to a broadening of the 
peaks. The peaks near the middle of the allowed band 
are wider than those near the edge, as clearly seen in 
Fig. 2c, which is due to the larger values of A oc dcu/dk'. 
Note that, although the integral over the absorption spec- 
trum @ monotonously increases with T m (a), the peak 
value of An(lu) has a maximum at T m (a) = Tm , i.e., 
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fm ) = 1/2. Analogous behavior of the 



FIG. 4: Transmission (a), reflection (b), and absorption (c) 
spectra of a superconducting photonic crystal calculated for 
three indicated values of a and the same set of parameters 
as in Fig. 2. Solid and dashed curves correspond to spectra 
calculated numerically and by using Eqs. (0. At a = 1 the 
absorption is zero, the corresponding line A(ui) = coincides 
with the abscissa and is indistinguishable. 



absorption peaks was observed for the two-dimensional 
photonic crystal consisting of spherical voids buried in 
a metal film [l2l |. If we take into account overlapping 
between peaks © with different m the maximum value 
of the absorbance can slightly exceed 1/2, as is seen in 
Fig. 2c. 



IV. CONCLUSION 

We have studied the photon dispersion law and opti- 
cal properties of a two-dimensional superconducting pho- 
tonic crystal formed by metallic cylinders arranged in a 
square lattice. The transmission, reflection and absorp- 
tion spectra are calculated as a function of the fraction 
a of superconducting electrons. In the spectral region 
within the allowed band the spectra exhibit oscillations 
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due to the interference effects. For the pure SC phase, 
a = 1 , the absorption vanishes, R{w) = l—T(u>), and the 
transmission peak values reach unity. With the transition 
from superconducting to normal metallic phase when a 
changes from 1 to the optical spectra exhibit remark- 
able modifications. The proposed approximate analytical 
description allows one to explain the key features of the 
dispersion and spectral modifications with changing a. 
The peaks in transmission and absorption spectra as well 



as the dips in the reflection spectra can be interpreted in 
terms of the resonance tunneling. 
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